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Homogenization of metamaterials is a crucial issue as it allows to describe their optical response 
in terms of effective wave parameters as e.g. propagation constants. In this paper we consider the 
possible homogenization of chiral metamaterials. We show that for meta-atoms of a certain size 
a critical density exists above which increasing coupling between neighboring meta-atoms prevails 
a reasonable homogenization. On the contrary, a dilution in excess will induce features reminis- 
cent to photonic crystals likewise prevailing a homogenization. Based on Bloch mode dispersion 
we introduce an analytical criterion for performing the homogenization and a tool to predict the 
homogenization limit. We show that strong coupling between meta-atoms of chiral metamaterials 
may prevent their homogenization at all. 

PACS numbers: 78.20.Ek, 81.05.Xj, 78.67.Pt, 42.25.Bs, 78.20.Bh 



I. INTRODUCTION 

At present chiral metamaterials (MMs) attract a great 
deal of interest of the scientific community since their po- 
tential to strongly affect the polarization state of light is 
promising for many applications. There are many mani- 
festations of this potential of chiral MMs as giant optical 
activity 1-3 , circular dichroismi~— and negative index of 
refraction^—, which are either considerably enhanced in 
comparison with the naturally occurring phenomena or 
cannot be observed in nature at all. 

By definition, the unit cell of a chiral material is not 
superimposable with its mirror image and any chiral ma- 
terial exhibits gyrotropy in general. However, a giant 
gyrotropy requires strong intrinsic resonances of the ele- 
mentary constitutive components, i.e. the chiral "meta- 
atoms" . For non-resonant structures the optical rotary 
power per meta-atom layer is usually vanishingly small. 
Even for resonant structures the achievable polarization 
rotation is often below the requirements of applications. 
This limitation can be overcome by stacking of several 
functional layers. Conventionally, light propagation in 
multilayered MMs is described by assigning effective wave 
parameters (EPs) to this Mtf ' u i 13 il , such as e.g. the 
propagation constant k or the effective refractive index 
introduced as n = kc/uj . Such EPs simplify the design 
and the description of functional MM devices. However, 
before assigning EPs it is of utmost importance to assure 
that the respective MM may be actually homogenized 1 -. 

When the mutual coupling between adjacent func- 
tional layers in a stack is strong, the convergence of e.g. 
the effective refractive index with the stack thickness to- 
wards its bulk value is poor—. Then, and this is a com- 
mon situation, the EPs depend on the MM slab thick- 
ness. Such effective properties have very limited sensed, 
since they cannot predict the response (for example, re- 
flectance and transmittance) of a MM slab with an arbi- 



trary thickness. The EPs are very rarely identical for a 
single and an infinite number of monolayers, though de- 
signs exist exhibiting this property 2 ^. Near field coupling 
has been investigated for the constitutive elements of 
meta-atoms 2 ^—, magnetic plasmon structure d 25 ' 26 and 
chiral MM bi-layers^. The latter work [27[ showed that 
the functional properties of the bi-layered MM with cou- 
pling between layers arc different from the properties of 
a monolayer. However, up to our knowledge, the gen- 
eral analysis of the coupling effects on the EPs of chiral 
metamaterials has not been conducted yet. 

The main purpose of this work is to provide a system- 
atic and comprehensive analysis of the influence of the 
coupling effects on the EPs of chiral MMs and on the 
chiral MMs homogenization. We introduce an analytical 
criterion of homogeneity based on the dispersion relation 
of Bloch modes. We show that several chiral MMs can 
be homogenized and in such case there is an optimal MM 
period, or in other words, an optimal distance between 
MM layers. We also show that in other cases MMs cannot 
be homogenized in the spectral resonance domain despite 
of a weakening of the coupling among adjacent layers by 
increasing their separation. 

We consider only effective wave parameters rather than 
any material EP s 20 ' 28 . Hence, effects caused by the meso- 
scopic structure of MMs (nonlocality, spatial dispersion, 
see e.g. j2Cl29j) are not in the focus of this contribution. 
The aim is rather to show under which conditions a mul- 
tilayered periodic chiral MM may be described as an ef- 
fectively homogeneous medium, similar to the investiga- 
tions performed for the fishnet M M 30 ' 31 and multilayered 
dielectric stack in^ 2 -. To avoid any complications with the 
tensorial properties of homogenized materials, the direc- 
tion of propagation is fixed along the axis pointing nor- 
mal to the MM surface and the responses of the media 
are probed upon illumination with left- and right-circular 
polarized waves only. Because of normal incidence the 
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eigenvalue problem reduces to a scalar problem. 

After having introduced the methodology of the re- 
search in Sec. II, we analyze as an example the twisted- 
cross MMi with respect to its homogenization in Sec. 
III. Additionally we introduce an analytical homogeniza- 
tion criterion. As a counterexample for a structure that 
cannot be homogenized, the twisted split-ring resonator 
MM22, is considered in Sec. IV. The importance of cou- 
pling effects for chiral MMs is emphasized in the Sec. V. 
Mechanical analogy of the periodic metamaterial is also 
provided there. The Conclusions section sums up the 
work. 



II. METHODOLOGY 

For an adequate description of a MM as homogeneous 
several requirements have to be fulfilled. 

i. The operating frequency must be below the first 
Bragg resonance 5i(fca) < 7T, where a is the meta- 
atom size in the direction of wave propagation. 
Otherwise, the MM constitutes a photonic crystal 
and its homogenization is meaningless^. 

ii. The lateral MM periodic must be such that 
the propagation of nonzero diffraction orders is 
suppressed 2 ^, i.e., the transverse size of a unit cell 
is required to be smaller than half the wavelength. 

iii. For each polarization state the light propaga- 
tion must be governed by a single (fundamental) 
Bloch mode only 2 - 9 having the smallest propagation 
losses. 

iv. The incident light should predominantly couple to 
this Bloch mod o 29 ' 30 . 

In general, situations may occur where light couples to 
a higher-order Bloch mode more efficiently than to the 
fundamental one. If this higher-order Bloch mode is more 
damped than the fundamental one, its impact becomes 
negligible above a certain thickness of the MM slab and 
EPs converge to those of the bulk. However, if the at- 
tenuation coefficients for all Bloch modes to which light 
can couple are comparable, then this leads to thickness- 
dependent EPs and results in the lack of homogeneity. 

As the tool to analyze such questions we rely on the dis- 
persion relation of Bloch modes. The respective disper- 
sion diagrams provide comprehensive information about 
the eigenmodes that are sustained by the MM. The oc- 
currence of a Bragg resonance can be extracted from the 
dispersion relation and occurs whenever 3?(fca) = 7r holds. 
From the analysis of the dispersion of the imaginary part 
of the wavenumber, ^s(k), we can evaluate the propaga- 
tion losses. So it is easy to infer about requirements (iii) 
and (iv) - number of dominating modes and their losses 
by considering in detail the dispersion relation. Thus, the 
Bloch modes dispersion diagram is an efficient mean of 




FIG. 1: (Color online). Metamaterial designs under consider- 
ation: TC (a) and TSRR (b). Wave propagation and meta- 
material stacking direction is along z-axis. 

MM homogeneity characterization. We would like to re- 
mark that in the case of chiral media with C4-symmetry 
there are two fundamental modes: right and left circular 
polarized (RCP and LCP). 

The chiral properties of a MM are determined through 
their effective refractive indices for the LCP and RCP 
waves the til and ur, respectively. The chirality param- 
eter is given by k = (tir — til)/2. Its real part 5R(k) 
describes the rotation of the polarization ellipse (or opti- 
cal activity )ii whereas the imaginary part S(k) governs 
the circular dichroism. 

The first homogeneity requirement implies that a 
smaller lattice period a in the propagation direction is 
beneficial. To check this statement we calculate the Bloch 
modes dispersion diagrams for both the twisted-cross 
metamaterials 1 (Fig.[T^) and twisted split-ring resonator 
metamaterials^ 3 . (Fig. [TL) . The lateral sizes of their unit 
cells are such that the homogeneity requirement (ii) is 
satisfied in the whole frequency range of interest. These 
examples were only chosen for the sake of proving the 
applicability of our approach and constitute referential 
unit cells that are currently under investigation. 

To simulate the optical response of MM slabs with 
various thicknesses we rely on the commercially avail- 
able CST Microwave Studio software 3 ^ that implements 
the finite-integral method in the frequency domain. The 
thin slab EPs are retrieved from the complex reflec- 
tion/transmission coefficients 15 , which we refer to as the 
standard method. The thick slab EPs, which are close 
to bulk parameters, are restored from the spatial depen- 
dence of the wave field with the wave propagation re- 
trieval method 3 - 5 -* 3 ^ . The Bloch dispersion diagram cal- 
culations are performed for infinite periodic media using 
a plane wave expansion method 3 ^ and taking into ac- 
count the dispersion of the metal permittivity. Details of 
the calculation method can be found elsewher o 29 i 30 ' 38 . 

Comparing the spectral simulations of finite structures 
performed with the CST Microwave Studio with those 
done with the plane wave expansion method we observed 
a steady blue shift in spectra of about 20 THz for the 
latter method. This shift that happens due to different 
numerical approaches is of minor importance for our con- 
clusions. 
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III. TWISTED CROSS METAMATERIAL 

The twisted-cross (TC) MMi (see Fig. QJi) is an ex- 
ample of a chiral MM that provides strong optical activ- 
ity. Its meta-atoms consist of two golden crosses embed- 
ded into a dielectric. One cross is 22.5 degrees twisted 
with respect to the other. Lateral unit cell sizes are 
a x = a y — 500 nm. The metallic crosses of thickness 
25 nm are separated with a dielectric spacer of 37.5 nm. 
The eigenfunctions of the TC are RCP and LCP because 
of the C4-symmetry of the MM. For our investigation we 
use the parameters of Ref. Q. However, in Ref. [l[ 
only a monolayer is characterized. To access the bulk 
properties of such a MM we stack multiple monolayers 
as presented in Fig. QJi. The meta-atom size in the prop- 
agation direction a is subject to variations. The crosses 
in each meta-atom are twisted by 45 degrees with respect 
to each other. Details can be seen in Fig. QJ, where the 
angles of consecutive layers are indicated. 

Note that, although the real unit cell size in propa- 
gation direction is 2a, the size a is actually the length 
scale of importance because it denotes the periodic dis- 
tance of consecutive functional layers. Since the optical 
response of the two parts of the unit cell are close to 
identical, we wish to regard this distance a as the period. 
In fact, calculating the dispersion diagrams for the struc- 
ture with the second double cross being identical to the 
first one, i.e. leading to a true period a, results qualita- 
tively in identical characteristics. This holds except for 
the case of dense stacking, where the material becomes 
achiral. The Bragg resonances for the present structure 
are caused by the reflection of light at the consecutive 
layers of the unit cell rather than by the boundary be- 
tween two unit cells. Bragg gaps due to the period 2a 
are not observed. Therefore, we will call the parameter 
a in the following the period of the structure. 

To challenge the third requirement of the homogene- 
ity we applied the standard retrieving (based on reflec- 
tion/transmission) to thin MM slabs and the wave prop- 
agation method to thick MM slabs. The aim of such a 
retrieval with different methods was to check the con- 
vergence of the effective parameters. We studied MMs 
with two different center-to-center separations, the dense 
(a = 125 nm) and the sparse stacking (a = 250 nm). The 
thin-slab case consisted of 1-3 meta-atom layers for both 
MMs. For the thick slab EPs retrieval, 48 monolayers 
were stacked, thus mimicking the bulk material. 

The structures exhibit completely different types of 
convergence with the slab thickness. For the dense stack- 
ing the convergence in the resonant region is very slow 
(Fig. [2^, and b). So we suspect that the TC is not homo- 
geneous there due to strong coupling between consecutive 
unit cells. The coupling can be reduced by separating the 
unit cells more. Indeed, for the sparse stacking, the ef- 
fective refractive index of the thin slab quickly converges 
to the bulk values that are retrieved with the wave prop- 
agation retrieval method (Fig. and d). This is a clear 
indication that such MMs might satisfy the homogeneity 
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FIG. 2: (Color online). Effective refractive index of left circu- 
lar polarized light of a twisted-cross metamaterial as a func- 
tion of frequency and for different numbers of layers that form 
the MM. Dense stacking (period 2a — 250 nm) $l(nLCp) (a), 
Sr(riLCp) (b); sparse stacking (period 2a = 500 nm) $t(n,LCp) 
(c), Sj(n_LCp) (d). Effective parameters are retrieved with the 
standard method 1 - (SM) for one (black line), two (green/grey 
line) and three (cyan/light grey line) monolayers and with the 
wave propagation retrieval method^ (WPRM) for 48 layers 
(red/dark grey dots). 

conditions. 

To check the limits where TC homogenization fails we 
calculated the Bloch modes dispersion diagram for in- 
creasing separation, as shown in Fig. [3] It is clearly seen 
that for the dense stacking (a = 125 nm) higher-order 
Bloch modes have a damping comparable with the fun- 
damental modes in the resonant region (Fig. [3£i). Slow 
convergence of the EPs as seen in Fig. b, confirms 
that light couples to the higher-order Bloch modes, thus 
homogeneity requirement (iv) is violated. 

For the sparse stacking case of a = 250 nm (Fig. [3)d) 
the continuum of the higher-order Bloch modes is sepa- 
rated from the fundamental ones in terms of losses. The 
higher-order modes are much stronger damped than the 
fundamental ones. Thus one may expect that only the 
fundamental modes govern the light propagation. More- 
over, the MM unit cells are still sub- wavelength (the res- 
onant response occurs at frequencies far below the low- 
est order Bragg resonance), so all four requirements are 
satisfied, and the TC MM can be considered as a homo- 
geneous. 

Increasing the period to a — 400 nm we observe a new 
peak in ^s(k) emerging at around 275 THz and which is 
separated from the lowest resonance by a local minimum 
in the losses at 250 THz (Fig. [St). The Bragg condi- 
tion 5ft(fcet) = 7T (the plot for 5R(fca) not shown here) is 
satisfied between 250 THz and 300 THz, therefore, this 
peak can be clearly identified as the first Bragg resonance. 
In the frequency region above the first Bragg resonance 
(^R(ka) > 7r) the properties of periodic arrangement of 
unit cells dominate and, according to requirement (i), 
homogenization fails. The frequency ranges where pho- 
tonic crystal properties prevail are shadowed in Fig. [3] 
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FIG. 3: (Color online). Bloch modes damping Sr(fc) dispersion 
spectra for twisted-cross MM for different separations: 125 
nm (a), 250 nm (b), 400 nm (c), 500 nm (d). Fundamental 
modes are shown as black line with circles and green/grey line 
with triangles, all the higher-order Bloch modes are shown 
in red/dark grey dots. The frequency region above the first 
Bragg resonance is shadowed. 



mental modes increase considerably for smaller separa- 
tion due to the increase of the metal filling fraction. 

The brief summary of the presented results is that 
there exists a range of optimal longitudinal periods to 
achieve the best homogeneity for the twisted-cross MM, 
which is limited by the influence of the higher-order 
Bloch modes from one side and the Bragg resonance 
from the other side. To determine it we plot the dif- 
ference AS(fca) = min i > 3 (Q ! (fc i a)) — max(Sj(fcia, k%a)) 
between the attenuation constants of the least damped 
higher-order Bloch mode and the most damped funda- 
mental mode (Fig. [4^,). We suggest to neglect higher- 
order Bloch modes in the consideration if the difference 
in damping per unit cell is A^(ka) > 0.5. In this case 
the intensity of the higher-order modes will be attenu- 
ated e times more over the period a, compared to the 
fundamental modes. In the resonant region between 200 
and 250 THz the homogeneity condition A^s(ka) > 0.5 is 
satisfied for a > 200 nm (Fig. |4^) . Outside the resonant 
region this condition is much more moderate and the TC 
can be homogenized for smaller values of a. 

The dependence of the minimum of the damping dif- 
ference min AS(fca) as a function of the period a is shown 
in Fig. Surprisingly, this dependence is not only 

monotonously rising but also close to linear. This fact 
points out the predictive potential of the used method- 
ology. Thus, it suffices to calculate the Bloch modes for 
only two periods in order to approximately predict the 
homogenization limit for the designed metamaterial. 

We also checked the homogenization criterion 
A3(fcet) > 0.5 for the split cube in cage negative index 
MM21. This MM shows very fast convergence of the EPs 
with the number of monolayers - the refractive index is 
the same from one monolayer to the infinite number of 
monolayers. In the whole frequency range of interest 
higher-order Bloch modes are well separated from the 
fundamental ones (AS(fca) > 1.5) that explains why the 
EPs of the split-cube in cage do not depend on the slab 
thickness. 



IV. TWISTED SPLIT-RING RESONATOR 
METAMATERIAL 



Increasing a even further to 500 nm we face the situa- 
tion where the first Bragg resonance overlaps with the in- 
trinsic resonance of the meta-atoms at 220 THz (Fig.[3Ji). 
Interestingly, the dip in the absorption spectrum at 210 
THz appears instead of the peaks for all other separa- 
tions (Fig. [3^, b and c). This resembles the phenomenon 
of electromagnetically induced transparenc y 39 ^ , when 
two coupled oscillators with a small detuning show an 
absorption dip at the frequency, where the individual os- 
cillators resonance is observed. 

If we compare the magnitudes of mode losses extracted 
from Fig. [31 we conclude that losses of the higher-order 
Bloch modes remain almost constant while decreasing the 
unit cell size. On the contrary, the losses of the funda- 



It is natural to make a hypothesis that for some meta- 
material designs with strong coupling between consecu- 
tive layers, the MM cannot be homogenized in the reso- 
nance region at all. In other words, A3(/ca) < 0.5 even 
for the large separation of meta-atoms where homoge- 
nization is prevented by the appearance of Bragg res- 
onances. To check this hypothesis we investigated the 
twisted split-ring resonator (TSRR) MM^ (see Fig.[lp). 
It consists of the golden U-shaped split-ring resonators 
embedded in a dielectric. The in-plane lattice constants 
are o>x 

= 885 nm. The split rings of thickness 60 
nm are separated by a dielectric spacer of 85 nm. All 
the geo metrical and material parameters are taken as in 
Ref. 33] except the period a. 
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FIG. 4: (Color online), (a) The difference of the propagation 
constants AQ(ka) — mirii> 3 [Q(kia)] — max[9(fcia, fea)] of the 
least damped higher-order Bloch mode and the most damped 
fundamental mode for various unit cell sizes a: 125 nm (black 
line), 150 nm (red line with squares), 200 nm (green line with 
circles), 250 nm (blue line with triangles up), 300 nm (or- 
ange line with triangles down), 400 nm (magenta line with 
stars), 500 nm (brown line with crosses), (b) The minimum 
of minA3 ! (fco) in (a) as a function of the period (red cir- 
cles). The straight line represents a linear fit of minA^(fca) 
on a. The dashed line corresponds to the homogenization 
limit A3(fa) = 0.5. 



We calculated the Bloch modes dispersion diagram 
for three center-to-center separations of the monolayers: 
a = 290 nm, a = 450 nm and a = 600 nm. One of the 
fundamental modes designated by the black solid line in 
Fig. [5] fails to satisfy the homogeneity requirements even 
for the sparsest stacking at a — 600 nm. Its losses are 
higher than those for a few of the higher-order modes. 
The other fundamental mode (Fig. [5] green line) exhibits 
a comparable 3(fca) with one of the higher-order Bloch 
mode for a = 290 nm and a — 450 nm. For the sparsest 
stacking a = 600 nm, this fundamental mode experi- 
ences a slightly lower, but still unsufficient for homoge- 
nization, damping, since A(3(fca)) = 0.09 < 0.5. A fur- 
ther increase of the lattice period to a = 600 nm shifts 
the real part of the wavevector into the first band gap 
(K(fca) = 7r) for frequencies above 80 THz, so the require- 
ment (i) is not met. We would hardly term a metama- 
terial homogeneous if only one eigenwave (for example, 
RCP) fulfills the homogenization requirements. We may 
conclude that the TSRR does not satisfy requirement (iii) 
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FIG. 5: (Color online). The loss spectrum for TSRR MM 
for different periods: 290 nm (a), 450 nm (b), 600 nm (c). 
Fundamental modes are shown as black lines with dots and 
green/grey lines with triangles, all higher-order Bloch modes 
are shown as red/dark grey dots. The frequency region above 
the first Bragg resonance is shadowed. 



around the resonance frequency for any period. This is 
caused by the strong coupling between its meta-atoms. 
Thus, we prove our hypothesis that there exist MMs that 
cannot be considered homogeneous at the resonance fre- 
quency region. 

It is important to mention that such MM keeps its op- 
tical functionality, i.e. affecting the state of polarization. 
However, one must not extrapolate the properties of a 
single or a few functional layers towards a stack of more 
functional layers. 



V. DISCUSSION 

The homogeneity criterion (i) requires a unit 
cell/ wavelength ratio to be as small as possible. How- 
ever, in this case, as shown above, in the resonance region 
requirements (iii) (single Bloch mode) and (iv) (efficient 
coupling of light to this specific mode) are frequently not 
met. When decreasing the unit cell size a, the fundamen- 
tal modes' damping increases (see Fig. [3]), so sooner or 
later it becomes comparable with the higher-order Bloch 
modes' damping. The higher-order Bloch modes then 
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FIG. 6: The sketch of the split-ring resonators stacking in 
the negative-permeability configuration (a) and in the chiral 
configuration (b). Arrows represent the direction of the wave 
propagation, electric and magnetic field vectors. Black lines 
in the lower pictures represent the cross-section view of the 
lines of the magnetic field created by SR A and penetrating 
SR B. 
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FIG. 7: Mechanical analogy of the metamaterial in the pho- 
tonic crystal regime (a), homogeneous (b) and inhomogeneous 
(c) medium. Small balls and springs represent host medium. 
Large balls represent resonant inclusions. 



have a comparable impact on the wave propagation in a 
MM slab as a direct consequence of the strong interaction 
of the meta-atom monolayers. 

This conclusion, which is basically applicable to any 
MM, is extremely important for chiral MMs where the 
coupling of monolayers is usually strong. For example, 
for a negative-permeability configuration of the split-ring 
(SR) resonator metamaterial in the side-coupled arrange- 
ment, where the wavevector lies in the SR plane, SRs arc 
weakly coupled due to the mutual inductance (magnetic 
field penetration) (see the sketch in Fig. |5^). The induced 
current in SR A creates a magnetic field that penetrates 
SR B. However, the number of the field lines penetrating 
SR B is not very large. On the contrary, for a chiral con- 
figuration, when the stacking direction is perpendicular 
to the SR plane (see Fig. [Sb) , the magnetic field of SR A 
strongly penetrates SR B. 

Most chiral metamaterials consist of split-ring res- 
onators or gammadions, where the meta-atoms are 
stacked similarly as in Fig. [5p. The other configuration, 
where the magnetic field of one meta-atom penetrates 
weakly into the other one, corresponds to the twisted- 
cross MM. It is basically an extension of the cut- wire 
pairs towards CVsymmetry. The currents excited in the 
double crosses create the magnetic field perpendicular to 
the stacking direction and are thus weakly penetrating 
into the adjacent monolayer (similar as in Fig. |6^i). 

Historically, e.g. by Newton, Bernoulli, Lagrange and 
others, wave phenomena were frequently illustrated by 
taking advantage of mechanical analogies. The historical 
overview of a periodic medium description with a spring- 
ball system can be found in Brillouin's book^i. We like- 
wise employ a mechanical spring-ball system to illustrate 
the coupling effect in our MM system, see Fig. [7J 

The host dielectric medium is represented by chains 
of small balls and springs. Meta-atoms are represented 
in Fig. [7J as massive oscillators having their own reso- 



nance frequency. When the distance between the meta- 
atoms is large, so that Sft(fca) > n, the periodic medium 
is in the photonic crystal regime (Fig. [7^,). The mas- 
sive oscillators periodically arranged in the system scat- 
ter the wave individually. Then constructive or destruc- 
tive interference leads to the formation of stop-bands and 
pass-bands. When the distance between meta-atoms de- 
creases, but coupling between monolayers is still weak, 
we get a kind of coupled resonator waveguide, where 
only nearest-neighbor coupling matters (Fig.[7)D). Such a 
waveguide has one fundamental mode (two in case of chi- 
ral metamaterials), which is dominating the light prop- 
agation. Thus the periodic medium may be homoge- 
nized, and the extraordinary properties of MMs originate 
from the meta-atom resonances. When the distance be- 
tween the meta-atoms is further reduced, the coupling 
gets stronger and the nearest-neighbor scenario ceases to 
hold. Now higher-order modes come into play (Fig. [7]:) 
and participate in the energy transfer together with the 
fundamental mode. So, this is the regime, when higher- 
order Bloch modes can propagate in the structure, and 
homogenization is not feasible. 

As can be inferred from the spring-ball model the effect 
of higher-order modes is negligible if there is no strong in- 
teraction between the meta-atoms. In this case the effec- 
tive properties are solely determined by the resonances of 
the individual meta-atoms. However, to ensure a decou- 
pling while the period has to remain small in comparison 
to the wavelength, i.e., avoiding photonic crystal effects, 
the constitutive elements of the meta-atoms have to be 
very small. Obviously, in this case the resonance features 
of the meta-atoms disappear for a fixed frequency. Such 
MM would hardly exhibit any extraordinary properties. 
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VI. CONCLUSIONS 

To sum up, it turns out that the possibility of homog- 
enization of chiral MM depends crucially on the coupling 
strength between meta-atoms. For weak coupling be- 
tween monolayers the MM can be homogenized. Then 
the optimal meta-atom separation can be identified by 
looking at the Bloch mode dispersion relation in the fre- 
quency domain below the Bragg resonance. The TC MM 
homogenization works best if the characteristic parame- 
ter a of the unit cell is as small as possible but the higher- 
order Bloch modes remain well separated from the fun- 
damental modes. In this contribution we propose the 
analytical homogenization criterion as AS(fca) > 0.5. In 
this respect it is helpful that the calculation of the dis- 
persion relation for two unit cell parameters a suffices to 
estimate the lower homogeneity limit of a metamaterial. 

For strong coupling between meta-atoms, the homoge- 



nization of MMs might be impossible at all. In this case 
the functionality of a MM, e.g. a stack of several mono- 
layers, must be characterized rather than retrieving EPs 
of one monolayer and extrapolating them to the proper- 
ties of a multilayered stack. 
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